Modeling the Q-cycle mechanism of transmembrane energy 

conversion 

Anatoly Yu. Smirnov 1 ' 2 and Franco Nori 1 ' 2 
1 Advanced Science Institute, RIKEN, 
Wako-shi, Saitama, 351-0198, Japan and 
2 Physics Department, The University of Michigan, Ann Arbor, MI 48109-1040, USA 

Abstract 

The Q-cycle mechanism plays an important role in the conversion of the redox energy into 
the energy of the proton electrochemical gradient across the biomembrane. The bifurcated elec- 
tron transfer reaction, which is built into this mechanism, recycles one electron, thus, allowing to 
translocate two protons per one electron moving to the high-potential redox chain. We study a 
kinetic model of the Q-cycle mechanism in an artificial system which mimics the bf complex of 
plants and cyanobacteria in the regime of ferredoxin-dependent cyclic electron flow. Using methods 
of condensed matter physics, we derive a set of master equations and describe a time sequence of 
electron and proton transfer reactions in the complex. We find energetic conditions when the bifur- 
cation of the electron pathways at the positive side of the membrane occurs naturally, without any 
additional gates. For reasonable parameter values, we show that this system is able to translocate 
more than 1.8 protons, on average, per one electron, with a thermodynamic efficiency of the order 
of 32% or higher. 

PACS numbers: 87.16.A-, 82.39. Jn, 87.16.D- 
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I. INTRODUCTION 



The energy produced by a biological system or by an artificial device must often be con- 
verted into a more stable form [l[ The thermodynamic efficiency and the quantum yield 
of this process is of prime importance for the overall performance of the energy transducer. 
This transducer consumes the energy of input particles, which move energetically downhill 
or just disappear in the process, and transfer this energy to another kind of particles moving 
energetically uphill. Here, the quantum yield {QY) is defined as the number of particles at 
the output of the energy transducer divided by the number of input particles. The efficiency 
of such device can be higher when the quantum yield is more than one, i.e., when a single 
input particle creates many output carriers. 



The 
photon 



generation of two or more electron-hole pairs 



excitons) by a single high-energy 
4|]. A similar situation takes place 



was observed in semiconductor nanocrystals 
in the bc\ complex embedded into the inner mitochondrial membrane as well as in the 
related complex bf, which mediates the electron transfer between the Photosystem II (PS 
II) and the Photosystem I (PS I) in the thylakoid membranes of plants and cyanobacteria jsj]. 
According to the generally accepted Q-cycle mechanism 6|-[8), the transfer of two electrons 
from a plastoquinol molecule PQH2 to plastocyanin (in bf complexes) is accompanied by an 
energetically-uphill translocation of four protons from the negative (N) to the positive (P) 
side of the membrane, resulting in a quantum yield QY = 2. 

We note that within the standard redox loop mechanism (see [s] and references therein), 
only two protons are transferred in parallel with the transfer of two electrons, implying that 
the quantum yield is equal to one. Hereafter, we primarily concentrate on the bf complex 



as a biological counterpart of our artificial system. Despite numerous studies 



10|4l3|, the 



physical mechanism of the Q-cycle in bc\ and bf complexes is not completely understood. 

In this work we analyze a simple model (see Fig. 1) mimicking the main features of the 
Q-cycle in the bf complex in the regime of ferredoxin-dependent cyclic electron flow |14j. I15|. 
In this regime, electrons cycle between the PS I and the bf complex, which are electronically 
connected by a pool of ferredoxin molecules (on the N-side of the membrane) and by a pool 
of plastocyanin molecules on the lumenal (P) side of the membrane. We treat these two 
pools as a source (S) and drain (D) electron reservoirs coupled to the electron-binding sites 
A and B, respectively. Besides the sites A and B, the membrane-embedded central complex 
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is comprised of sites L and H, which correspond to hemes b^ and bu of the complex bf. 
The sites A and H are assumed to be electronically decoupled as well as the sites B and L. 
A mobile shuttle Q (an analog of a plastoquinone molecule) diffuses inside the membrane, 
between the sites A and H (on the P-side) and the sites B and L (on the N-side). As its 
biological counterpart, the shuttle has two electron sites and two proton-binding sites. At 
its N-position, the shuttle takes one electron from the site A and another electron from the 
site H and transfers these electrons to the sites B and L. 

At the N-side, the shuttle also accepts up to two protons from the stromal (electrically- 
negative) proton reservoir and donate these protons to the lumenal proton reservoir at the 
P-side of the membrane. When the fully populated Q-molecule arrives at the P-side, one 
electron from the shuttle goes strictly energetically downhill, to the site B, whereas another 
one returns to the L-H chain to be loaded again on the shuttle. The origin of this bifurcated 
reaction [l6|, which occurs at the P-side catalytic center, remains unknown. 

Here, we explore physicochemical conditions wherein our artificial complex is able to 
translocate twice as protons as the number of electrons transferred energetically downhill, 
from the source S to the drain D. Protons are translocated energetically uphill, from the N- to 
P-side of the membrane. We aim at the explanation of the Q-cycle operation in this artificial 
complex. We examine a wide range of parameters allowing the efficient performance of the 
Q-cycle scheme. The functional principles of the Q-cycle in artificial systems can provide a 
better understanding of the Q-cycle mechanism in the natural bc\ and bf complexes. 

II. MODEL AND METHODS 
A. Components and states 

To simplify the problem, we divide the whole system of six electron and two proton- 
binding sites into four weakly-interacting subsystems: (i) the LH subsystem consisting of L 
and H sites (with energies El,£h)\ (h) the shuttle Q having two electron sites l e , 2 e (with 
energies E\ = e 2 = £q) and two proton sites l p ,2 p (with energies E\ = E 2 — Eq), (iii) the 
site A (with energy e^), and (iv) the site B (with energy eb) ■ 

The LH-subsystem is characterized by 4 microscopic states starting with the empty (vac- 
uum) state, \1)lh = \0l0h), and ending with the doubly-occupied system: \4)lh = |lz,l.ff)- 
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The electron and proton populations of the shuttle are described by 16 states, where 
|1)<2 = |0i e 02e0ip02p) is the vacuum state and |16}q = |li e l 2e lipl2p) is the state of the 
completely-loaded shuttle. Here, we use the notation a (1 Q ) for an empty (occupied) site 
a, where a = l e , 2 e for the electron sites on the shuttle, and a = l p , 2 p for the proton 
sites. The average populations of the sites A and B are denoted by {%) and (%}. Here 
(. . .) = ((^ol • • • |^o))t means double- averaging over an initial wave function \l/o and over 
a thermal distribution (. . .) of reservoirs and environment characterized by the common 
temperature T. 

We take into account strong Coulomb interactions between sites from the same subsystem. 
A Coulomb repulsion between electrons located on the sites L and H (LH-subsystem) is 
characterized by the energy ulh, whereas for the electrons and protons on the shuttle (Q- 
subsystem) we introduce the following parameters: U e (Coulomb repulsion between two 
electrons occupying the shuttle sites l e , 2 e ); U p (electrostatic repulsion between two protons 
located on the sites 1 P ,2 P ); and Uij = U ep (electron-proton Coulomb attraction on the 
shuttle). Here the indices i = 1,2 and j = 1,2 run over the electron and proton-binding 
sites of the shuttle, respectively. 

B. Electron and proton transitions 

Electrons can tunnel between the L and H sites as well as between sites belonging to dif- 
ferent subsystems. These tunnelings are described by relatively small tunneling amplitudes: 
A LH (L-H tunneling); A H q(x) (tunneling between the H site and the l e ,2 e sites; A L q(x) 
(for electron transitions between the l e , 2 e and L sites. The amplitudes Aaq{x) and Abq(x) 
describe the tunneling between the bridge sites A, B and the electron sites l e , 2 e on the 
shuttle. The amplitudes Aaq,A H q and A B q,A L q depend on the position x of the mobile 
molecule Q. The proton rates T^{x) (for transitions between the sites l p , 2 P and the N-side 
proton reservoir) and Tp(x) (for transitions between the proton sites on the shuttle and the 
P-reservoir) are also functions of x. 

The time evolution of x(t) is determined by an overdamped Langevin equation with terms 
describing walls of the membrane and a potential barrier, which impedes the charged shuttle 
to cross the intermembrane space [17|. The walls of the membrane are located at x = —Xq 
(N-side) and at X — Xq (P-side), so that the width of the membrane is equal to 2xq. 



4 



C. Equations and potentials 



Using the formalism outlined in Ref. 17H19j and based on Marcus transitions rates 20], 



we derive and numerically solve a set of master equations (see an Appendix) for an elec- 
tron distribution (Rm) over states of the LH-subsystem coupled to a system of equations 
for the probabilities (p M ) to find the Q-subsystem in the state These equations are 
complemented by equations for the populations (n^), (%) of the A and B-sites and by an 
overdamped Langevin equation for the shuttle position x(t). 

It is known [5] that charged ions, such as Mg 2+ and CI - , can easily cross a thylakoid 
membrane, which results in the equilibration of electrical potentials on both sides of a 
membrane. Therefore, the difference between the proton electrochemical potential /ip of the 
P-side, and the potential /jn of the N-side of the membrane 5] is mainly determined by the 
proton concentration gradient ApH : 

Hp- fi N ~ -ApH x (T/Tr) x 60 meV, (1) 

where T is the temperature of the reservoirs, and Tp = 298 K is the room temperature. 
It should be noted, however, that a surface potential Vs(x), which is positive on the N- 
side, Vs(— xq) = +Vn, and negative, Vs(xq) = — Vp , on the P-side of the membrane, was 



calculated for the bf complex of M. laminosus [21( with Vjv = 4.6 T and Vp = 5 AT (the 
Boltzmann constant kp = 1). This model assumes that there is a similar transmembrane 
potential (see Fig. 1), 

li <*> = - £ £f v, »- f £? l, » (2) 

with Vn — 120 meV and Vp ~ 140 meV, which correspond to the above-mentioned values 
of Vn and Vp at room temperature. All energies are measured in meV. In the presence of 
the surface potential, the energy levels of electrons and protons on the shuttle are shifted 
from their initial values Eqq and Eq depending on the shuttle's position x: 

£ q( x ) = e Qo ~ V s (x), 

E Q {x)=E Q0 + V s {x), (3) 

with e QN = e Q0 - V N , e QP = e Q0 + V P and E QN = E Q0 + Vn, E qp = E Q0 - V P . Correspond- 
ingly, the energy levels of the electron sites A and H, located near the N-side, are shifted 
down from their initial values: Ea = £ao — Vn, £h = £m ~ Vn] whereas the energies of the 
sites B and L, located near the P-side, are shifted up: Bp = £bo + Vp, El = £lo + Vp. 



III. "PASSENGER" SCENARIO OF THE Q-CYCLE 



Instead of searching over a multidimensional space of system parameters, we consider a 
reasonable sequence of events, which provides an optimal performance of the energy trans- 
ducer. In particular, we analyze a scenario where an electron transferred from the high- 
energy source reservoir S along the chain: S^-A^-Q^-B^-D to the drain D performs 
the main energetic function in the transfer of two protons from the N to the P-side of the 
membrane. Another electron traveling on the shuttle and recycled by the LH-system along 
the chain H — )■ Q — )■ L — > H plays a more passive role of a passenger, which is necessary to 
compensate a shuttle charge. 



According to the Marcus formula 



Kn' = |Ajj>|%/— — exp 
A,;,;' 1 



20] (see also an Appendix), 

{Si — Eii — A,;,;') 2 



(4) 



AX U ,T 

the rate Kui for an electron transition from the site i, with an energy £j, to the site i', with 
an energy Ey , has a maximum at £j = Ey + Xu> . Here Au> is the tunneling amplitude between 
the sites i and i', Xa> is the corresponding reorganization energy, which is due to electron 
coupling to an environment with temperature T. The shuttle can accept protons from the 
N-reservoir provided that the electrochemical potential of the N-side, /iat, is higher than the 
proton energy level on the shuttle. Protons move from the shuttle Q to the P-side reservoir 
if the energy of the Q-proton exceeds the P-side potential \ip. 



A. Sequence of events and energy relations 

We start with a situation when the empty shuttle (quinone) is near the N-side catalytic 
center, x = —Xq, and the Lif-system (analog of cytochrome b in the bf complex) is preloaded 
with one electron located presumably at the site H, which has a lower energy than the site 
L: El > Eh- The site A is also occupied with an electron taken from the electron source 
S. We have the following sequence of electron (e) and proton (p) transfer from and to the 
shuttle located near the N-side of the membrane: 

a) e : H Q, e h = e qn + X HQ . 

b) e : A -> Q, e a = e qn + U e + X A q- 

c) p : N -)> Q, jj, N > E QN - 2 U ep . 
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d) p : N ->Q, fi N > E QN -2U ep + U p . 

Here we have written relations between energy levels of electrons and protons which make 
possible these transfers. 

The shuttle loaded with two electrons and two protons travels to the P-side of the mem- 
brane, where the following sequence of electron and proton transitions occurs: 

e) e : Q ->■ B, e QP + U e - 2U ep = e B + X BQ . 

f) p: Q ->■ P, E QP - U ep + U P > fi P . 

g) p : Q -»■ P, P QP - £4 P > /i P . 

h) e : Q -> L, £ QP = e L + A LQ . 

Finally, an electron tunnels from the L to the H site: 

i) e : L ->■ #, e L = + A L h- 

The empty shuttle diffuses to the N-side of the membrane and the process repeats. We 
expect that two protons will be translocated from the N-side to the P-side of the membrane 
per one electron transferred from the source to the drain electron reservoir with a quantum 
yield QY = 2. 

Here we assume that, as in the case of the quinone molecule Q the shuttle populated 
with one electron (after step 1) does not bind a proton but accepts another electron (step 2). 
The doubly-reduced quinol is known to have a much stronger ability for binding two protons 
(see steps 3 and 4). At the P-side of the membrane the process presumably evolves in the 
opposite direction when the transfer of one electron from Q to the site B is accompanied by 
the unloading of two protons. In the absence of an attraction to two positive charges, the 
energy of the electron remaining on the shuttle goes up; thus, allowing its tunneling to the 
L-site. 

Here, an electron recycled by the LH-system plays a passive role of a shuttle's "passenger" 
since its transitions to and from the Q-molecule are not immediately accompanied by a 
proton transfer. Transitions of another electron, which is loaded to the shuttle from the 
source S (via the site A) and unloaded to the drain D (via the site B), are more closely 
coupled to the energetically-uphill proton translocation. 

It follows, from the relations a),h),i) in this section, that the recycling of one electron 
by the LH-chain, H — > Q — > L — > H, which lies at the heart of the Q-cycle, takes place if 
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the difference of surface potentials, 

AV = V s (-x ) - V s (x ) = V N + V P , (5) 

is of the order of the total reorganization energy along the recycling path: 

AV = X HQ + X LQ + X LH . (6) 

We see from the relations d) and g) in this section that the energetically-uphill proton 
transfer from N- to the P-side of the membrane is possible if the original energy of the 
proton on the shuttle, Eq , obeys the following inequality: 

fi N + 2 U ep -U P -V N > E Q0 > fi P + U ep + V P , (7) 

which can be true only for a sufficiently strong attraction potential, U ep , between electrons 
and protons on the shuttle 

U ep > up - (i N + AV + U p . (8) 

The relations a), b), e), h), i) in this section allow to estimate the original energies of the 
electron-binding sites counted, e.g,, from the level epo- 

£ao = £bo + 2U ep + Xaq + Xbq, 

£ho = £bo + 2U ep — U e + Xpq + Xhq, 

e Q0 — e B0 + 2U ep — U e + XpQ, 

£lo = £bo + 2U ep — U e + Xpq — Xlq. (9) 

We assume that the potentials of the electron source, /is, and electron drain, /x D , are of the 
order of the energies of the A and B sites, respectively: fis = £a, Hd = £b- Taking into 
account Eqs. ()9]) we obtain a relation for the source-drain energy drop, 

Us - HD > 2 U ep + Xaq + X BQ - AV (10) 

With Eqs. (16|HI) we obtain the following requirement for the energy difference between the 
source and drain electron reservoirs: 

/As - > 2 (lip - fjL N ) + 2U P + Atot, (11) 
where the combined reorganization energy, 

Atot = Xaq + Xbq + Xhq + Xlq + Xlh, (12) 
accumulates all losses along both electron transport chains: A — > Q — > B and H — > Q — > 
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B. Thermodynamic efficiency and quantum yield 



The thermodynamic efficiency r\ of proton translocation can be defined as 

t^m x Np (13) 

where Np is the number of protons translocated from the N- to the P-side of the membrane, 
and rip, is the number of electrons transferred from the source S to the electron drain D. 
The efficiency rj is proportional to the quantum yield 

QY = ^. (14) 
n D 

It follows from Eq. (II ip that, within the "passenger" scenario, the efficiency 77 of the electron- 
to-proton energy conversion can be estimated as 

77 = Vp-Vn x qy 

2(np-fi N ) + 2U p + \ tot V V ' 

This means that for a high electrochemical proton gradient, jip — hn S> U p + 0.5A tot , the 
efficiency i] has the maximum: rj = QY/2. Thus, in the ideal case, when QY = 2, the 
thermodynamic efficiency can reach the perfect mark: 77 = 1, when almost all electron 
energy is converted to the transmembrane proton-motive force. 



IV. RESULTS AND DISCUSSIONS 



A. Parameters 



In the model presented here an electron transport chain begins at the source reservoir S, 
corresponding to a pool of ferredoxin (Fd) molecules, which carry electrons from Photosys- 
tem I to the bf complex (see Fig. 1). The electron drain D is related to the high-potential 
chain of the bf complex comprised of the ion-sulfur protein (ISP), cytochrome /, and sol- 
uble plastocyanin (PC) molecules. Taking into account a redox potential of ferredoxin, 
E m = —0.41 V, and the fact that a redox potential of the ISP///PC chain, E m is in the 
range from 0.3 to 0.45 V Q, we estimate that the total energy drop between the source 
and the drain, Us — Hd, takes values from 710 meV up to 860 meV. 

We also assume that, as for the bf complex [211 ]. the surface voltage gradient, AV = 
Vjq + Vp = 260 meV, has been applied to the membrane, with a positive potential, Vg(— xq) = 
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Vn = 120 meV, at the N-side and a negative potential, V$(xo) = —Vp = —140 meV, at the 
P-side of the membrane. 

The system reaches its optimal performance when U ep = 610 meV. Hereafter, we assume 
that U e = U ep /2, U p = U ep /8, Ulh = 240 meV. We use the following values for the elec- 
tron transfer rates, 75 = 7d = 0.1 /ieV, proton transition rates, T N = T P = 2 /ieV, and 
for peak values of the electron tunneling amplitudes, Aaq(-xo) = Adq(xo) = 0.1 meV, 
Ahq(~ xo) = A^q(xo) = 0.06 meV. The coefficients 75 and 7d determine the rates of the 
electron transfer between site A and source S and between site B and drain D, respectively 
(see the definitions of 75 and 7^ in the Appendix). The rates T N and T P , describing proton 
transitions between the N-side reservoir and the shuttle as well as between the shuttle and 
the P-side of the membrane, respectively, are defined in the Appendix and in the subsection 
IE of the Supporting Information. For the diffusion coefficient, T> = T/(, we have the value, 
T> ~ 8 x 10~ 12 m 2 /s, which is close to experimental data for plastoquinone molecules in a 
lipid membrane [23[. Hereafter, we assume that the reorganization energies, corresponding 
to H-Q and L-Q transitions, are equal: Xhq = Xlq, and the same relation is true for the 
A-Q and B-Q transitions: Xaq = Xbq- 

To satisfy Eq. fl6]) we start with small values of the reorganization energies along the 
recycling pathway, Xlq = Xhq = 60 meV, and Xlh = 140 meV, which add up to the value 
of the surface gradient, AV = 260 meV. Here, we also have Xaq = 60 meV. For the driving 
force, fis — ^D = 850 meV, related to the electron energy drop in the bf complex, the system 
is able to translocate almost two protons, QY ~ 1.9, against the electrochemical gradient, 
Hp — = 150 meV, per each electron transferred from the source to the drain. This proton 
electrochemical gradient corresponds to the value ApH = —2.5 at room temperature, T = 
Tr = 298 K. The thermodynamic efficiency of the energetically-uphill proton translocation 
is about 33 % (77 ~ 0.33). 

B. Time evolution of a proton translocation process 

A proton translocation process is shown in Fig. 2, where we plot the time dependence 
of the total electron, uq = (rii) + (r^), an d proton, Nq = (Ni) + (N2), populations of 
the shuttle (Fig. 2b), together with the time-evolving position of the shuttle x(t) (Fig. 2a). 
Here, we also show (see Fig. 2c) the populations of the L-site, (til), and H-site, (nn), 
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as well as the average number of electrons (no) transferred from the source to the drain, 
as well as the average number of protons (Np) translocated from the N-side to the P- 
side proton reservoir. The brackets (. . .) are dropped in Figs. 2c and 2d for the notations 
of the populations, and throughout the paper, except in the Appendix. Data for Fig. 2 
are calculated at higher values of the reorganization energies: Xlq = \hq = 100 meV, 
\ LH = 250 meV, and X AQ = 100 meV. We assume that U ep = 610 meV, AV = 260 meV, 
Us — = 850 meV, and \ip — /iat = 150 meV. 

It can be seen from Fig. 2 that during 30 /xs the shuttle performs about 8 trips from the 
N-side (x = —x ) to the P-side (x = Xq = 2 nm) of the membrane and back, translocating in 
the process about 7 electrons (n D = 6.7) and 12 protons (Np = 12.2), with quantum yield 
QY ~ 1.8, and thermodynamic efficiency rj = 32%. 

At the N-side (see Fig. 2b) the shuttle accepts an electron from the initially populated 
site H and another electron from the source S (via site A) as well as two protons from the 
N-side proton reservoir. We note (see Fig. 2c) that site H is not completely depopulated. 
This means that there are events when both electrons occupying the shuttle arrive from the 
site A and the source S, shorting out the Q-cycle pathway. This leakage process increases 
the number of electrons transferred from the source to the drain (with the same number of 
protons), thus decreasing the quantum yield QY. 

Here we do not impose any additional restrictions, except a proper choice of energy levels, 
which are close to values given by Eqs. ( |7|9lfTTT) . with fig — 410, Ea = 465, Eh = 220, and 
Eqn = 160, Eqn = 982, at the N-side catalytic center (all energies are measured in meV). 
For the P-side center we use the following energies: Eqp = 420, Eqp = 722, for electrons 
and protons on the shuttle, Ep = —495, \ip> = —440, for the high-potential redox chain, 
and el = 360 meV for the recycling pathway. This choice of energy levels makes the H-to- 
Q electron transition (at the N-side) much easier than the A-to-Q electron transfer, since 
£p — Eqn < Xhq, whereas Ea — £qn 3> Xaq- Moreover, the S-to-A electron transition is also 
hampered since the energy level of the A-site, Ea, is higher than the potential of the source, 
PS- 

On arrival at the P-side of the membrane the shuttle donates an electron to the B-site 
and, finally, to the drain D. Two protons move to the P-side proton reservoir (Figs. 2b and 
2d). It is evident from Fig. 2c that another electron from the shuttle Q goes to the L-site 
(see the small spike at the bottom of Fig. 2c). This electron is rapidly transferred to the 
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H-site, and the empty shuttle returns to the N-side. 

Figures 2b, 2c, 2d illustrate the bifurcated reaction which occurs at the P-side. Here, one 
electron from the shuttle Q goes to the high-potential (and low-energy) chain, Q — > B — > 
D, while another electron (a passenger) returns to the LH-system for recycling (along the 
pathway Q — )■ L — >■ H). No additional gate mechanisms are required for this reaction. An 
escape of the first electron from Q to B, followed by the transition of two protons to the 
P-reservoir, increases the energy of the remaining electron to the level Eqp = 420 meV, 
which is of order of the L-site energy, Ep = 360 meV, but is much higher than the energy of 
the B-site, Bp = —495 meV. Furthermore, the site B is probably occupied with an electron 
taken from the drain since jip > £b- These two factors strongly suppress the leakage of the 
second electron from the shuttle to the high-potential chain. 

C. Effects of the proton electrochemical gradient 

In Figure 3 we show the numbers of protons, Np, and electrons, np>, transferred across 
the membrane, as well as the quantum yield, QY, and the power-conversion efficiency, 
rj, as functions of the proton electrochemical gradient, ftp — fiN (measured in meV). The 
graphs are plotted for two sets of reorganization energies: (i) Xlq = Xhq = 100 meV, 
Xlh = 250 meV, Xaq = Xpq = 100 meV (blue curves) (ii) Xlq = Xhq = 200 meV, 
Xlh = 400 meV, Xaq = Xbq = 200 meV (green dashed curves). Other parameters, such as 
U ep = 610, us — \id = 850, AV = 260 (in meV), are the same as in Fig. 2. 

We numerically calculate the output of the system (Np, np,, etc.) at the end of the 
stochastic trajectory x(t) (with the duration t = 100 lis) and average results over 10 tra- 
jectories. For each value of the shuttle's position x(t) we solve a set of master equations 
(1AlllA5|Alip . which have been averaged over electron and proton reservoirs as well as over 
fluctuations of the environment coupled to the electronic degrees of freedom. This can be 
done since the electron and proton transitions are much faster than the mechanical motion 
of the shuttle. 

It follows from Fig. 3 that more protons, Np ~ 45, and electrons, rip, ~ 25, are transferred 
across the membrane at lower proton gradients, \ip — /ijy — 150 meV with a higher quantum 
yield, QY > 1.8. However, the thermodynamic efficiency is higher, rj ~ 39%, at larger 
proton gradients, fijy — \ip — 200 meV, where Np ~ 40 and QY ~ 1.6. These numbers 
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are for the set (i), with smaller values of the reorganization energies (see blue continuous 
curves in Fig. 3). A stronger electron-environment interaction, described by the set (ii) of 
reorganization energies, significantly reduces an energetically-uphill proton flow with almost 
no impact on the electron current (see green dashed curves in Fig. 3). In this case the 
quantum yield drops to almost one, which means that the recycling pathway (via the LH- 
system) is practically closed. 

D. Effects of the surface potential gradient 

Figure 4 demonstrates the performance of the system as a function of the surface potential 
gradient, AV = Vn + Vp, at a fixed difference between the N-side and P-side potentials, 
Vp — Vv = 20 meV. Here we choose two sets of system parameters. The first set (see the blue 
continuous curves in Fig. 4), with U ep = 610, fis — fiD = 850, and fip — fiN = 150 meV, was 
considered before. The second set (green dashed curves in Fig. 4) is characterized by a higher 
electron-proton attraction potential, U ep = 800 meV, and a higher source-drain difference, 
fi-S — fiD = 1220 meV. At these parameters the system can translocate protons against the 
electrochemical difference fip — fiN — 300 meV, which is related to the proton concentration 
gradient ApH = —5 at room temperature. The reorganization energies correspond to the 
set (i) described before. 

It follows from Fig. 4 that the system performs very well in a wide range of surface voltage 
gradients, from AV ~ 230 meV up to AV ~ 350 meV, translocating more than 40 protons 
and 20 electrons (in 100 fis). These numbers are almost the same for both sets of parameters 
(see Figs. 4a and 4b). The quantum yield QY monotonically goes down, from QY ~ 2 at 
AV ~ 200 meV to QY ~ 1.5 at AV ~ 350 meV, no matter what set is chosen. The 
power-conversion efficiency rj also decreases with increasing AV. However, r] is higher for 
the second set of parameters. For example, at AV = 260 meV, the efficiency is about 31% 
for the first set, and rj ~ 47% for the second set (with QY = 1.9) where both electron and 
proton electrochemical gradients are higher. This trend is consistent with Eq. (fl5|) derived 
for the passenger scenario. 
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V. CONCLUSIONS 



We have theoretically examined a model of the protonmotive force generation by the 
Q-cycle mechanism mimicking the operation of the bf complex in the thylakoid membranes 
of plants and cyanobacteria. We concentrate on a simulation of the regime of ferredoxin- 
dependent cyclic electron flow, where the bf complex translocates up to two protons (a 
quantum yield QY = 2) across the membrane per one electron transferred from the electron 
source (ferredoxin pool) to the electron drain (pool of plastocyanin molecules). This model 
includes two electron and two proton sites on the shuttle Q (an analog of a plastoquinone 
molecule), diffusing inside the membrane, as well as two electron sites, A and B, connecting 
the electron sites on the shuttle to the source and drain reservoirs. The recycling of an 
electron by the cytochrome 6, which forms the basis of the Q-cycle, can be described by 
adding two electron sites L and H corresponding to hemes and bu of the bf complex. 

We have derived and numerically solved a set of master equations for the populations of 
the electron and proton-binding sites together with a Langevin equation for the position of 
the shuttle. Within a reasonable scenario and in the presence of the surface potential, we 
have determined the conditions which are necessary for the efficient translocation of protons 
across the membrane. We have found that the system is able to transfer, on average, about 
1.8 protons per one electron (QY = 1.8) with a thermodynamic efficiency of the order of 
32% against the transmembrane proton gradient fip — fijy = 150 meV at the source-drain 
difference of electron potentials /is — /id = 850 meV. These values of the electron and proton 
gradients are closely related to experimental values for the bf complex. No conformational 
gating is necessary for the bifurcation of the electron transfer reaction at the P-side catalytic 
center, where one electron goes to the drain and another electron returns to the LH-chain, 
to be loaded on the shuttle again. We have studied the performance of the model as a 
function of the proton electrochemical gradient and the surface potential. It is shown that 
the system demonstrates even better results, with a quantum yield of the order of 1.9 and an 
efficiency of the order of 47%, when both the source-drain difference and the proton gradient 
are higher. 
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Appendix A: Master equations for electron-driven proton transfer across a mem- 
brane 



Here we briefly outline a set of master equations describing the process of electron-driven 
proton translocation across a membrane. As we mentioned before, the total system is 
characterized by the average populations of the A and B-sites, (ua), (ub), as well as by 
four states of the LH-subsystem with electron distributions (Rm) (M = 1,...,4) and 16 
states of the Q-subsystem (electrons and protons on the shuttle) with distributions (p M ) 

(/!=1,...,16). 

The time evolution of the LH-system is governed by the equation 



(R 



M 



N 



+ J2^m h n(Rn) 

N 



(Al) 



with the following relaxation matrix 



Imn — 7mtv + Imn + 1m n- 



(A2) 



Here the rate 



Imn - \Alh\ 



71 



x { \(M\a[a H \N)\ 2 + \{N\a\a H \M)\ 2 } 



x exp 



(A3) 



describes the L-to-H electron transitions. Hereafter, a' a , a a refer to the creation/annihilation 
operators for an electron on the site a. For protons on the site the creation/annihilation 
operators are denoted by A^,Ap. The rate jff N (and a similar rate ■y M %) is related to the 
electron transfer between the L (or H) sites and the sites l e , 2 e on the shuttle, 



Imn 



7T 



\(M\a L \N)\ 2 exp 
\(N\a L \M}\ 2 exp 



\<J*\ai + a 2 \v}\ 2 x 



(Umn — + Xlq) 2 1 



(Q>MN + + ^Lq) 2 

4A^T 



(A4) 
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where 01,02 are operators of the electron-binding sites on the shuttle, VLmn = Em — E^ 
are frequencies of the LH-system (h — 1, ks — 1), and u^ u is the frequency spectrum of the 
coupled electron-proton states on the shuttle Q. 

For the distributions of the 16 states of the Q-system we derive the following equation 



where 



(/v> = - i% (p») + E i% 



1 /if 1 /if 1 ' /if 1 1 /if 1 1 /if 1 f /if 1 1 fj,u ' 



(AS) 



(A6) 



Components of this relaxation matrix can be written as 



7^ Q = |A^ Q | 



7T 



|(/x|ai + a 2 k)| 2 exp 
+ + a 2|/i)| 2 exp 
with a similar matrix 7^7, and 



4A AQ T 



Aaq^ 
(1 - n A 



X 



(A7) 



ry L Q 
I /if 



\{v\a-y + a 2 \fi}\ 2 (Rn) exp 
+ |(/i|ai + a 2 \v)\ 2 (R M ) exp 



Alq|%/^^|(M|o l |7V)| 2 
V Xl Q T mn 

{^liv + &MN + Alq) 2 



X 



/if 



4A LQ T 



4A LQ T 



(A8) 



with a matrix 7^, which is similar to 7^. The proton transitions to and from the shuttle 
are described by the rate 



7 " Q = rVK/il^ + A 2 \u)\ 2 [l - F N {u^)\ + 



and by a similar rate 7^. Here 



T N \{v\A 1 + A 2 \n)\ 2 F N 



(, xp I — y: — I - ! 



(A9) 



(A10) 



is the Fermi distribution of the protons in a— reservoir (a = N,P). 
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The average population of the A-site is governed by the equation 



IA 



AQ | 



(n A ) = Is [fs{e A ) - (n A )] + 
\(n\ai + a 2 |^)| 2 x 



7T 



[1 -n A )exp 
-{n A )exp 



to 



+ € A + X 



AQj 



\21 



CO 



+ e A - X 



AQ) 



\21 



4A AQ T 



(All) 



where a coefficient 75 (or 7^) describes the electron transitions from the A (or B) site 
to the source S (drain D) electron reservoir characterized by a Fermi distribution with 
electrochemical potentials /is or /ic (a = S,D), 

-1 



/«(e) 



cxp 



T 



+ 1 



(A12) 



A similar equation takes place for the population (%). 

We solve the rate equations ( 1A1IIA5IIA11I) for both the distributions (Rm), {p^) and for 
the populations (n A ) and (ub), together with an overdamped Langevin equation for the 
mechanical motion of the shuttle, 



dU w 
dx 



dU, 



((n 1 + n 2 -N 1 -N 2 r)^+^ 

dx 



(A13) 



where ( is the drag coefficient, £ is a Gaussian fluctuation source with zero mean value, 
(£) = 0, and with a correlator (£(£)£(£')) = 2(T5(t — t'). The potential U w (x) confines the 
shuttle between the membrane walls, and the potential U c h(x) prevents the charged molecule 
Q from crossing the lipid core of the membrane (for details see Ref. 17| ). We note that the 
tunneling amplitudes A A q, Ahq and the proton rate depend on the distance between the 
shuttle (with a coordinate x) and the N-side catalytic center located at x = — x , whereas 
the tunneling amplitudes Abq, Alq and the proton rate Tp depend on the distance between 
the shuttle and the P-side catalytic center located at x = x . 
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FIG. 1: Simplified schematic of the Q-cycle mechanism in the regime of cyclic electron flow. At 
the N-side of the membrane the Photosystem I (PS I) provides the source reservoir S with high- 
energy electrons. Via the bridge site A, the source S delivers electrons to the shuttle Q, which also 
accepts electrons from the site H and protons from the N-side proton reservoir. At the P-side of 
the membrane the shuttle Q gives away electrons to the site L and to the drain reservoir D (via 
the bridge site B). In this process, two protons move to the P-side proton reservoir. From the site 
L electrons return to the site H to be loaded later on the shuttle Q. The drain reservoir D transfers 
low-energy electrons back to the Photosystem I. The surface potential Vs, which is positive at the 
N-side and negative at the P-side of the membrane, is shown here with circled plus and minus 
signs. 
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FIG. 2: (a) Time evolution of the position x{t) of the shuttle diffusing between the membrane 
walls located at x = ±2 nm; (b)-(d) a sequence of electron and proton transfer reactions at the 
following set of parameters: fis — [J>d = 850, fip — hn = 150, U ep = 610 meV, and AV = 260 meV. 
The total proton population Nq of the shuttle [blue continuous line in (b)] almost coincides with 
the electron population uq marked by the dashed green line. The bifurcated electron transfer 
reaction takes place at the P-side of the membrane (at x = 2 nm) where one electron moves from 
Q to the site B and the drain D [see, e.g., a step down for uq in (b) and a step up for np> in (d) 
at the moment t ~ 3 /is]. At almost the same time another electron moves to the site L [see a 
barely visible dashed red line in (c)] and rapidly proceeds to the site H [blue continuous spike in 
(c)]. Two protons are unloaded from the shuttle to the P-side of the membrane as follows from the 
step for Np [blue continuous curve in (d)] . It can be seen from (d) that the number of protons Np 
translocated to the P-side (blue continuous line) is nearly twice as large as the number of electrons 
np, moved to the drain (dashed green line). 
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FIG. 3: (a) Number of protons Np transferred to the P-side of the membrane; (b) number of elec- 
trons rip, moved to the drain; (c) their ratio (quantum yield QY); (d) the thermodynamic efficiency 
rj as functions of the transmembrane proton gradient Lip — fiN (measured in meV) for different cou- 
plings to the environment (in meV): weaker couplings (i) Xlq = 100; Xlh = 250; Xaq = 100 (blue 
continuous curves); and stronger couplings (ii) Xlq = 200; Xlh = 400; Xaq = 200 (green dashed 
curves). Other parameters are the same as in Fig. 2. The data are averaged over 10 stochastic 
realizations. Each realization lasts for 100 fis. These graphs demonstrate the ability of the system 
to translocate protons against the gradient up to 200 meV with an efficiency r/ up to 40%. 



21 



60| r 




150 200 250 300 350 

Surface potential gradient AV (meV) 



FIG. 4: Dependence of the output indicators of the system, the number of protons translocated to 
the P-side, Np (a), the number of electrons moved to the drain, no (b), the quantum yield QY(c) 
and the efficiency r] (d), on the surface potential gradient AV = Vn+Vp across the membrane. Two 
kinds of curves correspond to two sets of energy values: (1) ^5 — fj,p> = 850; \xp — \xn = 150; U ep = 
610 meV (blue continuous curves), and (2) fig — \id = 1220; [ip — ^N = 300; U ep = 800 meV (green 
dashed curves) and for Xlq = 100, Xlh = 250, Xaq = 100 meV. It can be seen from this figure 
that the system is operational in the range of the surface potentials from AV = 230 meV up to 
AV = 350 meV. In agreement with Eq. (|15p the complex works more efficiently (with r/ ~ 50% and 
QY ~ 2) for the second set of parameters where both electron and proton gradients are higher. 
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